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ABSTRACT 

^ : 

Oh' We search for ongoing major dry-mergers in a well selected sample of local Bright- 

Q I est Cluster Galaxies (BCGs) from the C4 cluster catalogue. 18 out of 515 early-type 

\^ • BCGs with redshift between 0.03 and 0.12 are found to be in major dry-mergers, 

' which are selected as pairs (or triples) with r-band magnitude difference Sm-[ < 1.5 

Ci . and projected separation Tp < 30 kpc, and showing signatures of interaction in 

the form of significant asymmetry in residual images. We find that the fraction of 

^ I BCGs in major dry-mergers increases with the richness of the clusters, consistent 

^ , with the fact that richer clusters usually have more massive (or luminous) BCGs. 

' We estimate that present-day early-type BCGs may have experienced on average 

I ^ 0.6 (inicrge/0.3 Gyr)^^ major dry-mergers and through this process increases their 

■ luminosity (mass) by 15% (tniorgc/0.3 Gyr)"-^ (/niass/0.5) on average since z — 0.7, 
O^l ' where imorgo is the merging timescale and /mass is the mean mass fraction of compan- 

. ion galaxies added to the central ones. We also find that major dry-mergers do not 

' seem to elevate radio activities in BCGs. Our study shows that major dry-mergers 

I involving BCGs in clusters of galaxies are not rare in the local Universe, and they are 

. an important channel for the formation and evolution of BCGs. 



Key words: galaxies: elliptical and lenticular, cD - galaxies: clusters: general - galax- 
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1 INTRODUCTION 

The Brightest Cluster Galaxies (BCGs) are at the most lu- 
minous and massive end of galaxy population. They are usu- 
ally located at or close to the centres of dense clusters of 
galaxies based on the X-ray observations or gravitational 
lensing observations (e.g., Jones & Forman 1984; Smith et 
al. 2005). Most of them are dominated by old stars with- 
out prominent ongoing star formation. They are much more 
luminous than other galaxies in the clusters, and have dis- 
tinct surface brightness profiles from other galaxies (see Vale 
& Ostriker 2008 and references therein). And yet the lumi- 
nosities of BCGs appear to have rather small scatters. As a 
result, they have been proposed as a distance indicator (e.g. 
Postman & Lauer 1995). Because of the unusual properties 
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of BCGs, their formation and evolution has long attracted 
particular attention. 

Galactic cannibalism as a way to build up the BCGs 
was proposed more than three decades ago (e.g., Ostriker & 
Tremaine 1975; White 1976; Vale & Ostriker 2008). How- 
ever, as pointed out by several authors (e.g., Merritt 1985; 
Lauer 1988; Tremaine 1990), cannibalism of satellite galaxies 
may not be enough to account for the luminosity growth of 
massive BCGs due to the small luminosity ingestion rate. If 
this is correct, then other mechanisms, such as major merg- 
ers may play an important role (Lin & Mohr 2004; Vale & 
Ostriker 2008). Since most central galaxies in clusters have 
little cold gas, the mergers are not expected to trigger ma- 
jor episodes of star formation and thus they are likely dry 
mergers. 

Recent studies from both numerical simulations and ob- 
servations indicate that a large part of stellar mass in lu- 
minous early-type galaxies form before redshift z ~ 1 — 2, 
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and later dry mergers play an important role in their stel- 
lar mass assembly (Gao et al. 2004; De Lucia & Blaizot 
2007; van Dokkum et al. 1999; Tran et al. 2005). Ruszkowski 
& Springel (2009) investigated the effect of dry-mergers on 
the scaling relations of BCGs in simulations. Observational 
studies on the properties of luminous early-type galaxies 
provide evidence that they have experienced dry-mergers 
(e.g., Lauer et al. 2007; Bernardi et al. 2007; von der Linden 
et al. 2007; Liu et al. 2008, hereafter Paper I). Theoretical 
investigations suggest that half the mass of a typical BCG 
may be assembled via accretion of smaller galaxies, i.e., by 
minor mergers at z < 0.5 (e.g., Gao et al. 2004; De Lu- 
cia & Blaizot 2007). Recently, Bernardi (2009) studied their 
evolutions in the size- and velocity dispersion-stellar mass 
correlations since z < 0.3, and concluded that early-type 
BCGs grow from many dry minor mergers (see Discussion). 
However, Whiley et al. (2008) using the ESQ Distant Clus- 
ter Survey, found little evolution of BCGs since z ~ 1, at 
least within a metric circular aperture of 37 kpc. A consensus 
about the evolution of BCGs is yet to emerge. 

Recent modelling suggests that major dry-mergers play 
a significant role in forming luminous, intermediate-mass, 
early-type galaxies (e.g. Naab et al. 2006) at redshift z ^ 0.5. 
Whether the same scenario applies to the much more lumi- 
nous BCGs, and more generally how the occurrence of ma- 
jor dry-mergers depends on the environments, is still some- 
what controversial. Merritt (1985) argued that massive ma- 
jor mergers are more likely to occur in large groups than in 
massive clusters, since the merger rate is a steeply declining 
function of the velocity dispersions of a virilized system. 

However, recent observational studies indicate major 
dry-mergers do occur in groups and clusters of galaxies. For 
example, Mulchaey et al. (2006) and Jeltema et al. (2007) 
reported some examples of dry-mergers involving central 
galaxies in intermediate-redshift groups. Tran et al. (2008) 
reported an observational analysis of supergroup SG 1120- 
1202 at z ~ 0.37, which is expected to merge and form a 
cluster with mass comparable to Coma. They argued that 
the group environment is critical for the process of major 
dry-mergers. Rines et al. (2007) reported a very massive 
cluster (CL0958-I-4702) at moderate redshift {z = 0.39), in 
which a major dry-merger is ongoing to build up a BCG. 
Using the Sloan Digital Sky Survey (SDSS) data, Mcintosh 
et al. (2008) identified 38 major mergers of red galaxies from 
845 groups of galaxies at z ^ 0.12 identified using the halo- 
based group finder of Yang et al. (2005), and showed that 
centres of massive groups are the preferred environment for 
major dry-mergers and mass assembly of massive early-type 
galaxies. 

All these studies show that major dry-mergers are im- 
portant in the formation of BCGs at moderate redshifts. Our 
own visual inspections of 249 merging pairs identified from 
luminous early-type galaxies (Wen, Liu & Han 2009) also 
hint that some of these mergers involve BCGs. The purpose 
of this work is to perform a search for major dry-mergers 
in a well-selected local early-type BCG sample with a quan- 
titative method and then estimate their major dry-merger 
and luminosity (mass) increase rates. We also investigate 
the dependence of the merger fraction as a function of the 
richness of the clusters. 

The structure of the paper is as follows. In §2 and §3 
we describe our sample and identifications of major dry- 



mergers involving BCGs. We present our main results in §4 
and finish with a summary and discussion in §5. Through- 
out this paper we adopt a cosmology with a matter density 
parameter f2m — 0.3, a cosmological constant flA ~ 0.7, 
and a Hubble constant of Ho ~ 72 kms^^Mpc"^, i.e., 
h = //o/(100kms-^Mpc-^) = 0.72. 



2 THE SAMPLE 

The C4 cluster catalogue (Miller et al. 2005) was constructed 
from the SDSS spectroscopic data in the parameter space of 
position, redshift, and color. The C4 algorithm identifies the 
BCG from the spectroscopic catalogue within 0.5h~^ Mpc 
of the central galaxy at the peak of the C4 density field 
("mean" galaxy). The early (DR2) version of the C4 cat- 
alogue incorporates the SDSS photometric catalogue to se- 
lect the brightest cluster galaxy within lh~^ Mpc in order to 
avoid missing BCGs in spectroscopic data due to fiber col- 
lisions. The catalogue is relatively complete in the redshift 
range 0.03 ^ z ^ 0.12, including 98% of X-ray-identified 
clusters and 90% of Abell clusters in the same region (Miller 
et al. 2005) . We thus limit our data within this range of red- 
shift, obtaining 643 clusters after we reject 16 duplicates 
(Bernardi et al. 2007). 

However, as pointed out by several previous works (e.g., 
Bernardi et al. 2007; von der Linden et al. 2007), some bright 
stars have been mis-classified as BCGs, and about 1/4 of 
BCGs show late-type features (e.g., distinct spiral arms, 
dust features). We thus re-checked each cluster and iden- 
tified the brightest galaxy within lh~^ Mpc of the cluster 
center as the BCG. We then select a sample of early-type 
BCGs from this catalogue by excluding artifact stars and 
late-type galaxies. Furthermore we require the BCGs should 
be in clusters with richness larger than 10, and have colour 
typical of early-type galaxies, g — r > 0.7 (Strateva et al. 
2001) in the SDSS model magnitudes. In the end we obtain 
a sample of 515 early-type BCGs in the local Universe. Fig- 
ure [T] shows the distributions of redshifts of the clusters and 
g — r colours of BCGs, which has taken into account the 
extinction by the SDSS, and k-correction to z = 0.1 (Blan- 
ton & Roweis 2007). The cluster redshifts range from 0.03 
to 0.12, with a maximum around 0.08; the g — r colours of 
BCGs peak around 0.9 with a small RMS scatter of around 
0.06, consistent with their stellar populations being quite 
old and uniform. 



3 IDENTIFICATIONS OF MAJOR 

DRY-MERGERS INVOLVING BCGS 

Dry-mergers between gas-poor early-type galaxies are gen- 
erally more difficult to identify than those mergers involv- 
ing late- type galaxies (Naab et al. 2006). The latter often 
develop prominent tidal tails dotted with star-forming re- 
gions. However, dry-mergers have more subtle features, such 
as broad "fans" due to ejected stars, asymmetries in in- 
ner isophotes and/or sometimes diffuse tails (Rix & White 
1989; Combes et al. 1995). While major dry-mergers usually 

^ The richness is defined as the number of member galaxies within 
lh~^ Mpc of the central galaxy. 
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Figure 1. The top panel shows the redshift distribution of red- 
shift Zci for the 515 C4 clusters while the bottom panel shows the 
colour g — r distribution of their early-type BCGs. 



have more prominent morphological signatures than minor 
dry- mergers (Le Fevre et al. 2000), they still exhibit much 
weaker merging features than wet-mergers. 

We adopt the following criteria to select candidates 
of ongoing major dry-mergers involving the BCG in our 
sample. First, galaxies with projected distance from central 
BCG must satisfy rp < 30 kpc. We search for close pairs both 
automatically and visually. In the automated search, a lower 
limit of the projected distance, 7 kpc, is imposed, which cor- 
responds to a separation of 3" at z — 0.12. This restriction 
is necessary because the photometry becomes unreliable for 
galaxies with smaller angular separation (Mcintosh et al. 
2008; Wen, Liu & Han 2009). Pairs or triples with separa- 
tions smaller than 7 kpc are searched visually. Second, they 
should have magnitude diflerence with their corresponding 
central BCGs < 2 magnitudes in the SDSS r— band model 
magnitude. The magnitude difference applied here is slightly 
larger than the value (1.5) that will be eventually used to 
select major-mergers. The reason is that the SDSS magni- 
tudes for BCGs are somewhat inaccurate due to the problem 
of sky background subtraction (Paper I). Third, both galax- 
ies in a pair must have g — r > 0.7 to ensure the merging 
galaxies are early-type (Strateva et al. 2001). In total, there 
are 49 apparent close pair candidates at 7 < rp < 30 kpc 
satisfying these three conditions (pairs with rp < 7 kpc are 
discussed below). 

Most previous studies impose an additional requirement 
of a small difference in the line-of-sight velocities 5v^ (e.g.. 




Figure 2. The r— band residual magnitudes mros versus the 
asymmetry factor A. The black points are the pairs in Wen et 
al. (2009, Fig. 5 in their paper). The 36 merger candidates with 
7 < rp < 30 kpc are shown with red circles. Four additional merg- 
ers with rp < 7 kpc are shown with red triangles. The horizontal 
and vertical dotted lines illustrate the criteria of Wen et al. e.g., 
rrires < 20 mag and A. > 0.45 used to select major mergers. The 
histograms for A and m^es for the Wen et al. (2009) sample are 
shown at the right and bottom for major-mergers with strong 
interactions (solid histogram) and those classified as weak inter- 
actions or due to chance alignment (dashed histogram). 



Svr < 500 km s^^) besides a small projected distance rp to 
avoid chance alignment. However, SDSS spectroscopic sur- 
vey is severely incomplete on very small angular scales due to 
fiber collision (Mcintosh et al. 2008; Wen, Liu & Han 2009). 
Hence it is not possible to apply this criterion to the SDSS 
data. An alternative method is needed to identify merging 
systems from the projected pairs and triples. 

Although, as mentioned above, the interaction features 
(e.g., broad plumes at the outskirts, short tidal tails, bridges 
or asymmetries in inner isophotes) are weak in mergers 
of early- type galaxies, they can nevertheless be seen more 
clearly from the residual image after we subtract a smooth 
and symmetric model for each galaxy in an image (Bell et 
al. 2006; Mcintosh et al. 2008; Wen, Liu & Han 2009). Such 
a procedure has been described in detail by Wen, Liu & 
Han (2009). For completeness, we repeat the essentials of 
the method below. 

First, we extract the SDSS r— band images of all 49 
close pairs and triples. The sky backgrounds are subtracted 
precisely from the corrected frames with the method of Pa- 
per I. Briefly, the background is subtracted as follows; all 
the objects including target sources in the corrected frame 
were flrstly masked to obtain a background-only image. 
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Second-order Legendre polynomials are then used to fit the 
smoothed background-only image to obtain an accurate sky 
background model. After background subtraction, the GAL- 
FIT package (Peng et al. 2002) is then used to construct a 
smooth symmetric Sersic (Sersic 1968) model for every early- 
type galaxy in the projected pair or triple. Stars and fainter 
galaxies with distances closer than 27?9o from the centre of 
BCGs (Here i?9o is the radius containing 90% of Petrosian 
flux from the SDSS catalogue) are modelled with a double- 
Gaussian point-spread-function (Stoughton et al. 2002) and 
Sersic function respectively. Objects outside 2i?9o in the ex- 
tracted images are masked. All models are convolved with 
the point-spread-function; we then obtain the best model 
by minimizing the with the sky-subtracted image of un- 
masked pixels. 

A fitted magnitude for each target galaxy and a resid- 
ual image for each pair are obtained afterwards. Notice that 
the fitted magnitude is superior because it can separate the 
flux in the overlapping region of pairs and correct the over- 
estimate in the sky background in the SDSS pipeline (Paper 
I). Wen, Liu & Han (2009) improved the method of Conselice 
et al. (2000) to calculate the asymmetry factor for residual 
images of galaxy pairs. Essentially it measures the difference 
between pixels and those symmetric pixels with respect to 
the centre and major axis of each galaxy (see eqs. 8-10 in 
their paper). The pixels in overlapping regions are treated 
separately (see their Figure 4) . The asymmetry factor ^ ~ 
is for a galaxy pair without any interaction feature ideally. A 
large A means a stronger asymmetric interaction. In order 
to find an efficient criterion, they compared the asymme- 
try factors of two test samples of residual images selected 
through careful visual inspections: a sample with distinct 
interaction features and another almost without. Both are 
not contaminated by other objects within 2i?9o of target 
galaxies. The test showed that a close pair can be classi- 
fied efficiently as a merging system with distinct interaction 
features when the residual images have an r— band magni- 
tude Tros < 20 mag and an asymmetry factor A > 0.45. 
In this work, we performed further extensive tests and find 
that both the asymmetry parameter A > 0.5 and residual 
magnitude rres < 19.5 can select strongly interacting merg- 
ers efficiently (see the two histograms in Figure [2J . These 
criteria are only slightly different from those used by Wen, 
Liu & Han (2009). The criteria of rros < 19.5 mag is almost 
the same as requiring the fiux ratio of the residual image to 
the initial image to be larger than 2% (/rcs//ini > 2%, see 
the last column in Table [!}. 

If we limit major mergers as having a fitted magnitude 
difference Srrir < 1.5, we find that 36 out of 49 close pairs 
or triples satisfy 7 < rp < 30 kpc. Furthermore, 14 out of 
36 have residuals rros < 19.5 mag and asymmetry factors 
A > 0.5 (see Figure[2]). Others may be weak interacting sys- 
tems or pairs in chance alignment; two examples with small 
residuals are shown in Figure|3] We added 4 mergers (C1035, 
C1176, C3150, C3157) with rp < 7 kpc (e.g., mergers with 
close double nuclei) identified visually - these systems have 
fiber magnitude difference Snii < 1.0, and can be classified 
as merging systems if we apply the same criteria as merg- 
ers at 7 < rp < 30 kpc. The r— band observational images 
and corresponding residual images for these 18 mergers in- 
volving BCGs are shown in Figure |3l Their corresponding 
parameters are listed in Table [1] 
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Figure 4. Two typical example galaxies which after subtraction 
show small residuals and small asymmetry factors. They are likely 
pairs produced by either weakly interacting with each other or 
chance alignment. Each image is 80" X 80". 

Bernardi et al. (2003) derived the luminosity function 
of early-type galaxies. In the r-band, the magnitude for an 
L* galaxy (derived from a Gaussian fit) is A'l, = —21.09 
(taking h = 0.72 instead of their h = 0.70) . The BCGs 
in our sample have luminosities ranging from L« to roughly 
^ 10Z/«, with a mean luminosity of ~ 5Z/«. For C4 2089, the 
merging pair has luminosities 11. 2L, and 5.6Lt. Similarly, 
for C4 1055, the triple merging galaxies have luminosities 
7.91/, , 6.61/, and 0.461/, ; the two most luminous components 
have identical radial velocities, and thus there is little doubt 
that they are physically associated. Mergers appear to occur 
even among the most luminous galaxies in clusters (see H4.2p . 
If all the stars coalesce to form a single giant galaxy, then 
the luminosities will be as high as 16.81/, and 15L,; such 
galaxies are rarely seen. This suggests that a substantial 
amount of stars of the companion are tidally stripped to 
form an "envelope" around the central galaxy, a point we 
return to in H4.ll 



4 RESULTS 

4.1 Merger rate and stellar assembly of BCGs 

We have identified 18 (~3.5% of the total) major dry merg- 
ing pairs (or triples) involving the BCG in 515 C4 clusters 
at 0.03 ^ z 0.12 . It has usually been assumed that physi- 
cally bound pairs will merge due to dynamical friction. The 
merger timescale is, however, uncertain since it depends on 
a number of factors, including the separation, relative veloc- 
ity and mass ratios of galaxies. Kitzbichler & White (2008) 
noticed that previous studies on merger rates based on pair 
statistics have yielded a wide variety of results (see their 
paper for a review). This diversity can be attributed to dif- 
ferences in the pair definition and in the timescales adopted. 
They thus used the Millennium Simulation to calibrate the 
merging timescale as a function of the stellar mass and pro- 
jected separation. They find an average merging timescale 
(appropriate for our stellar mass and the Hubble constant): 

Wsc = 2.2 Gyr ( ^ + , (1) 

^ ^ 50kpc V5.6 X IQiOMey/ V 8/ ^' 

where Af, is the total stellar mass, and z is the (cluster) red- 
shift. We take a mass-to-light ratio of M*/Lr « 5Mq/Lq 
(Rines et al. 2007; Patton et al. 2000), yielding a mean total 
stellar mass of M, ~ 8.5 x 10^^ Mq for our massive merg- 
ers. We estimate the merger timescale for each merger by 
this formula. We find the merger timescale of our mergers 
ranges from 0.04 to 0.55 Gyr with a mean value of 0.3 Gyr 
and a standard deviation of 0.17 Gyr (see Column 11 in Ta- 
ble[T]). For comparison, if we apply the formula of Masjedi et 
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Figure 3. The colour images and corresponding residual images in the r— band for the 18 identified major mergers. Each image is 
80" X 80". The SDSS-C4 name for the DR2 version is marked at the top left corner of each colour image. 



al. (2006) to calculate the dynamical timescale of our merg- 
ers, we obtain a similar timescale of ~ 0.3 Gyr. However, a 
somewhat shorter dynamical friction timescale of ~ 0.1 Gyr 
is obtained if we apply the formula of Patton et al. (2000) , 
although the dynamical friction formula may be less appli- 
cable to major mergers (Binney & Tremaine 1987). We also 
note that Bell et al. (2006) obtained a shorter timescale of 
~ 0.15 Gyr from simulations for massive early-type merg- 
ers. It seems that the average merging timescale is likely 
uncertain by a factor of ~ 2. 

To evaluate the number of mergers BCGs has experi- 
enced since (for example) z = 0.7, we need to understand 
how the dry- merger fraction evolves as a function of redshift. 
Parameterizing the evolution of pair fraction in the form of 
(l-|-z)™, several previous studies (e.g., Le Fevre et al. 2000; 
Patton et al. 2002; Conselice et al. 2003; Lin et al. 2004; 
De Propris et al. 2005; Kartaltepe et al. 2007; Lotz et al. 
2008) derive a positive slope in the range of m ~ — 3 for 
all mergers (involving dry or wet mergers). Lin et al. (2008), 
for the first time, estimated the merger fraction of red galax- 
ies (the majority of them should be dry-mergers) decreases 
with redshift slightly with m — —0.92 ± 0.59. This suggests 
that dry-mergers may preferentially occur in the local Uni- 
verse. In contrast, a recent work by Khochfar & Silk (2008) 
suggests a nearly constant dry-merger rate at 2 ^ 1. 

If we adopt an average merger timescale, and assume 
a constant fraction (18/515 ~ 3.5%) of BCGs are in ma- 
jor dry-mergers over the 6.3 Gyr interval since redshift 



0.7, then a present-day early-type BCG will have experi- 
enced 6.3/fmorgc X 18/515 ~ 0.7 {t 

merge /0.3 Gyr) major 
dry-mergers. On the other hand, if the fraction evolves as 
(1 -f 2)-of2±o.59 (Ljjj 2008), then the number of merg- 

ers changes to ~ 0.5(tmerge/0.3 Gyr)~^. Below we adopt an 
average number of mergers of these two evolution scenar- 
ios, 0.6 (tmcrgc/0.3 Gyr)~^. In any case, the uncertainty in 
the merger timescale is much larger than that due to the 
evolution of the merger rate. 



Our merger pairs/triples have an average magnitude 
difference of ~ 0.7, which corresponds to a ~ 1 : 2 lumi- 
nosity ratio. In the merging process, a substantial fraction 
of the companion galaxy may be tidally stripped, and so 
not all the mass of the companion galaxy will be added to 
the main one (Yang et al. 2009). We assume a fraction of 
/mass of the companion galaxy is accreted to the primary 
galaxy. Thus in each merger, the central galaxy increases 
its luminosity (or mass) by 25%(/mass/0.5). It follows that 
a present-day BCG has on average increased its luminos- 
ity by 15% (tnicrgc/0.3 Gyr)~^(/niass/0.5) from z = 0.7 at a 
rate of 2.5% {t 

merge /0.3Gyr)-i(/mass/0.5) per Gyr. Thus a 
non-negligible fraction of the stellar mass of a present-day 
early-type BCG is assembled via major dry-mergers since 
z = 0.7. 
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Table 1. Basic parameters for the 18 identified major dry-mergers involving BCGs 



C4 ID Mtot R.A.(J2000) Dec.(J2000) Zsp 5v^ m,. Mr tmcrgo A m,.cs /rcs//i, 

(mag) (kpc) (kms^^) (mag) (mag) ( Gyr) (mag) (%) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) 



1000 0.0864 -23.36 

1004 0.0814 -23.30 

1011 0.0903 -23.91 

1026 0.0908 -24.16 

1035* 0.0444 -22.95 

1055 0.0836 -24.03 

1060 0.1170 -23.29 

1176* 0.0737 -23.90 

1304 0.0986 -22.69 

1364 0.0708 -22.75 

2049 0.0751 -23.89 

2089 0.0652 -24.15 

2179 0.0938 -22.24 

3059 0.0406 -22.47 

3150* 0.0456 -22.05 

3157* 0.0635 -24.10 

3166 0.1112 -24.04 

3167 0.0915 -23.36 
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0.0866 

0.0863 15.94 90 
0.0814 

13.18 

0.0903 

22.92 
31.50 

0.0908 17.78 
0.0446 

0.0442 1.64 120 
0.0836 

0.0836 16.99 
8.93 

0.1170 8.13 
0.0738 

0.0736 6.37 60 
0.0987 

0.0986 20.72 30 
0.0710 

0.0707 21.27 90 

0.0751 7.44 
0.0652 

16.95 

0.0938 17.17 

0.0406 16.97 
0.0456 

2.21 

0.0632 

0.0638 4.97 180 
0.1112 

29.63 

0.0915 

0.0919 27.66 120 
0.0911 14.55 90 



15.42 -22.62 0.32 

15.44 -22.60 

14.90 -22.99 0.27 

16.09 -21.80 

14.87 -23.27 0.39 

15.51 -22.63 
16.38 -21.76 

14.50 -23.65 0.28 

15.06 -23.09 

16.11" - 0.04 
16.16" 

14.61 -23.34 0.28 

14.81 -23.14 

17.70 -20.25 

15.87 -22.91 0.17 
16.81 -21.97 
17.06" - 0.11 
17.62" 

16.12 -22.23 0.50 

16.84 -21.51 

15.06 -22.49 0.50 

16.49 -21.06 

14.05 -23.64 0.13 

15.52 -22.17 

13.65 -23.71 0.27 

14.40 -22.96 

16.71 -21.53 0.47 
16.80 -21.44 

14.26 -22.00 0.43 

14.94 -21.32 

16.43" - 0.06 
17.09" 

16.37" - 0.08 
16.89" 

14.93 -23.72 0.49 
16.11 -22.54 

15.88 -22.29 0.55 

15.94 -22.23 
16.25 -21.92 



0.75 17.98 5.2 

0.78 19.04 2.5 

2.71 17.79 5.1 

1.15 18.36 2.3 

0.69 18.46 2.1 

0.62 18.39 2.1 

3.49 19.09 4.2 

1.53 17.34 3.3 

0.87 18.55 7.4 

0.82 18.54 4.9 

0.77 18.04 2.2 

0.86 17.92 5.1 

0.57 18.43 10.6 

1.12 17.62 4.4 

0.71 17.72 4.8 

1.23 18.03 2.3 

0.53 19.29 3.3 

1.92 19.01 3.5 



Note: Col:(l). The C4 cluster ID (DR2) for the Clusters. The four visually-selected ones (indicated by *) have Vp < 7kpc (see JSj. 
Col: (2). Redshifts of identified merging pairs (or triples). Col: (3). Total absolute magnitude of mergers in the SDSS r— band. Col: (4). 
R.A. (J2000) of the component galaxy in a merger. Col:(5). Dec. (J2000) of the component galaxy in a merger. Col:(6). Spectroscopic 
redshift of the component galaxy in a merger. Col:(7). The projected distance, rp. Col:(8). The line-of-sight velocity difference, Svr, if 
available. Col:(9). Extinction-corrected r— band fitted apparent magnitude of component galaxy in a merger. For four mergers with 
rp < 7 kpc (marked with 'a'), we compare their fiber magnitudes. Col:(10). r— band absolute magnitude of component galaxy in a 
merger, corrected for extinction (by SDSS) and the k-correction (using the KCQRRECT algorithm of Blanton & Roweis 2007). Col:(ll). 
The merger timescale estimated by the formula of Kitzbichler & White (2008). Col:(12). The asymmetry factor for the residual image. 
Col:(13). r— band apparent magnitude of the residual image. Col:(14). The fiux ratio of residual and initial image within 3Re (Here Re 

is the fitted effective radius from our model). 



4.2 Dependence on environments 

In the hierarchical structure formation model, clusters of 
galaxies form through mergers of sub-clusters and small 
groups; in the process central galaxies are expected to grow 
together with clusters. This scenario suggests that the evo- 
lution of BCGs may depend on the environments. We thus 
examine the relation between the fraction of BCGs involved 
in major dry-mergers and the richness of the parent clusters. 



Figure [5] shows the results. The fraction of merging 
BCGs appears to increase as the clusters become richer. A 
simple linear fit in the form of fraction = A log richness -I- B 
yields A = 0.11 ± 0.04, B = -0.12 ± 0.05 (shown as the 
solid line), and so the significance is not very high due to 
the large error bars. The higher frequency of more massive 
dry-mergers in richer clusters will lead to more luminous (or 
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Figure 5. The top panel shows the major dry- merger fraction 
as a function of the richness of the parent clusters. Poisson er- 
rors are shown. The straight line shows the best regression fit 
fraction = (0.11 ± 0.04) log richness - (0.12 ± 0.05). The bottom 
panel shows the BCG luminosities versus the richness for 85 BCGs 
(open circles) in Paper I and 18 merging BCGs presented in this 
work (solid circles). 



more massive) BCGs. This is empirically seen in observa- 
tions, as illustrated in the bottom panel of Figure [5] 



5 SUMMARY AND DISCUSSION 

In this work, we have searched for BCGs involved in ma- 
jor dry-mergers in a sample of early-type BCGs selected 
carefully from C4 cluster catalogue (Miller et al. 2005). 18 
(~3.5%) major dry-mergers involving BCGs have been iden- 
tified out of 515 clusters at 0.03 ^ z ^ 0.12. These major 
mergers are selected as pairs (or triples) with r-band magni- 
tude difference Sirir < 1.5 and projected separation rp < 30 



kpc, and showing signatures of interaction in the form of 
significant asymmetry in residual images. 

To double check whether our results are affected by the 
choice of the cluster (C4) catalogue (based on DR2 of SDSS), 
we have also searched for major dry-mergers involving BCGs 
in the sample of 625 clusters refined by von der Linden et 
al. (2007) from the unpublished DR3 version of C4 clus- 
ters. Their sample selects preferentially early-type BCGs 
and discards clusters with very few member galaxies. For 
this sample, we find 27 (27/625~4.3%) BCGs are involved 
in major dry-mergers, entirely consistent with the fraction 
(3.5%) found here for the C4 cluster catalogue. It should 
be emphasized that the goal of our method is to identify 
merging pairs with strong interaction features and calculate 
their true faction in BCGs. They should be real, physically 
bound pairs. However, some physical bound pairs may not 
yet show strong interaction signatures. Our estimated frac- 
tion should thus be considered as the lower limit of physical 
pair fraction in BCGs. 

From the identified major-merger candidates, we con- 
clude that a present-day BCG may have experienced 
0.6 [t 

merge 

/0.3Gyr) ^ major dry-mergers and increased 
their mass by 15% (tmcrgo/0.3 Gyr)~^(/inaaa/0.5) on average 
from z = 0.7. This fraction is comparable to the mass in- 
crease predicted by mostly minor mergers in semi-analytical 
galaxy formation models of De Lucia & Blaizot (2007) , if all 
mass of companion galaxies are added to BCGs (/mass ~ 1) 
as De Lucia & Blaizot (2007) assumed. However, the mass 
increase we found here is due to major-mergers rather than 
minor majors. Bernardi (2009) concluded that the evolu- 
tion of BCGs in the velocity dispersion-stellar mass relation 
favors the growth of BCGs via many dry minor-mergers. 
However, the evidence is not very strong (see her Fig. 8). A 
more detailed comparison between theory and observations 
is desirable. 

The major uncertainties in our estimates arise due to 
the calibration of the merging timescale and the fraction of 
the satellite galaxy that will be added to the central host. 
For the former, we used the calibration of Kitzbichler & 
White (2008) derived from the Millennium simulation. How- 
ever, it remains to be seen whether their calibration applies 
well to the dense environments such as clusters of galax- 
ies. The parameter /mass is also somewhat uncertain. The 
stellar mass not accreted to the central host may form the 
envelopes in cD galaxies and/or the intracluster light (Lin & 
Mohr 2004; Rines et al. 2007; Paper I and reference therein). 
We notice, however, that if all the mass is assembled into 
the final BCGs, some of them may be too bright; the re- 
sulting luminosity function at the bright end may be incon- 
sistent with the observations which show little evolution at 
L > 2.5L* (e.g.. Brown et al. 2007; Cool et al. 2008). We 
conclude that a substantial fraction of the companion mass 
must be stripped, consistent with the presence of envelopes 
in giant cD galaxies (e.g.. Paper I). 

We have also checked the radio properties of BCGs to 
investigate whether major mergers trigger the phenomenon 
of active galactic nuclei. We cross-correlate our sample with 
the radio data from the FIRST survey (Becker, White & 
Helfand 1995). There are 94 (~18.3%) sources with radio 
emission out of 515 BCGs; these fractions are consistent 
with the statistics of Best et al. (2007) and Croft, de Vries 
& Becker (2007). In comparison, there are 3 radio sources 
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(~16.7%) out of 18 merging BCGs, i.e., there is no difference 
on radio properties between major-merger and non-major- 
merger BCGs within the statistical uncertainties. 

Mcintosh et al. (2008) identified 38 major mergers of 
red galaxies from 845 groups at z ^ 0.12, out of which 
18 mergers (~ 2.1%) are between a central galaxy and a 
satellite galaxy. Six (ClOOO, C1004, ClOll, C1060, C1304, 
C3157) of our 18 mergers are also in their sample of 38 ma- 
jor mergers; all six are in the top ten most massive groups 
among their major mergers. In other words, we are probing 
the most massive end of their SDSS "groups" in terms of 
mass range from 10^^-^ Mq to IO^'^Mq. They find that more 
massive mergers tend to appear in richer systems, a conclu- 
sion confirmed by our studies. In addition, they find that 
the merger frequency increases as a function of halo mass, 
which may explain why their merger fraction (~ 2.1%) is 
slightly lower than ours (3.5%). The mass increase rate we 
find for the BCGs is roughly consistent with their results 
(2%-9% per Gyr under similar assumptions). We conclude 
that major dry-mergers in BCGs are not rare events in the 
local Universe and they are an important way to assemble 
luminous BCGs. 
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